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Abstract 

The spin-boson model (or the dissipative two-state system) is a model for the study 
of dissipation and decoherence in quantum mechanics. The spin-boson model with 
Ohmic dissipation is an integrable theory, related to several other integrable systems 
including the anisotropic Kondo and resonant level models. Here we consider the prob¬ 
lem of computing the overlaps between two ground states corresponding to different 
values of parameters of the Ohmic spin-boson Hamiltonian. We argue that this can 
be understood as a part of the problem of quantizing the mKdV/sine-Gordon inte¬ 
grable hierarchy. The main objective of this work is to analyze how the Anderson 
orthogonality affects the Yang-Baxter integrable structure underlying the theory. 
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1 Introduction 

The occurrence of infrared (IR) divergences is a central issue for quantum field theories and 
condensed matter systems which possess gapless excitations BE]. Boundary Conformal Fields 
Theories (CFT) in two-dimensional space-time provide an opportunity to gain useful insights 
into the problem. In the simplest set up, with the Euclidean geometry of the half-inhnite 
plane (x, y) in which x < 0 is treated as a space coordinate, the vacuum states corresponding 
to different conformal Boundary Conditions (BC), say “1” and “2”, are orthogonal. More 
precisely, for a large but finite system of space size L their overlap tends to zero as a power 
^ defining an orthogonality exponent ^ 21 - This vacuum overlap can be interpreted as a one- 
point function of the BC changing operator 021 [3], and therefore its contribution to the spectral 
sum of the two-point function { 0 \i{y) 02 i{^)) vanishes, as well as individual contributions of 
all overlaps involving conformal descendant states. However, the combined contribution of the 
states organized in conformal towers turns out to be hnite for the inhnitely large system and 
gives rise to the scale-invariant two-point function (C>2i(2/)Oi2(0)) = IA 21 P The latter 

coincides with the ratio Z 2 i{y)/Zii, with ^21 standing for a partition function of the half- 
inhnite system with BC “1” everywhere except on a part of the boundary of length y, where 
BC “2” is imposed, whereas the denominator Zn is a partition function of the system with 
BC “1” is imposed along the whole boundary. There is generally a normalization ambiguity 
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Figure 1; An universal part of the overlap modulus |(f22|f^i)| can be extracted from the 
partition function of a classical statistical system with an inhomogeneous boundary. In this 
picture fii and fi 2 two numerically different sets of values of boundary couplings. 


of the numerical coefficient A 21 , and usually it can be set as 1. The situation becomes more 
interesting for the so-called boundary flows, i.e., for a class of two-dimensional quantum held 
theories in which conformal invariance is broken only by BC depending upon a set of couplings 
= {fii, fi 2 ...). Contrary to the conformally invariant case, the partition function of the system 
schematically visualized in Fig|T] is a complicated function of Euclidean time y. However, its 
y —)■ -l-cxo asymptotic is expected to have the form 

^^ = \A2i\^ + , (1.1) 

where AE '21 = Eq{^ 2 ) - is the difference in the ground state energies for the different 

sets of boundary couplings. The orthogonality exponent ^21 and the prefactor A 21 are functions 
of Hi and (-12 which virtually dehne an universal (i.e. independent on details of both IR and 
ultraviolet (UV) regularizations) part of the vacuum overlap 

{H 2 I ) oc A 21 as L ^ 00 . (1.2) 

In common nomenclature, the modulus of this overlap is referred as to the (ground state) hdelity. 
Throughout this paper, with some abuse of conventional terminology, this term will be used to 
denote the scaling function A 21 as well as its generalization. The generalization deals with the 
vacuum-vacuum matrix elements of a bare (unrenormalized) boundary held 0{y) characterized 
by a certain scaling exponent D{0), so that fll.2p is substituted by 

(H 2 I 0(0) I Hi) oc A2 i(0) L-d 2 i{o) ^ ^^ 3 ^ 

dehning both the IR exponent d2i{0) and the hdelity A 21 (O) (here e —)■ 0 is the lattice spacing, 
i.e., the UV regulator). 

The signihcance of the study of hdelities is that it may help to better understand univer¬ 
sal aspects of the dynamics after a local quantum quench in quantum impurity models [OHS]. 
Such models are used to mimic the behavior of small interacting quantum mechanical systems 
coupled to an external environment. In some cases, they display universality which can be 
described in terms of the boundary hows (see ref. [1] for review of applications of the bonndary 
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flows in condensed matter physics). Here we will discuss the so-called spin-boson model (or the 
dissipative two-state system) which is a paradigm model for study of dissipation and decoher¬ 
ence in quantum mechanics [5] . In the case of Ohmic dissipation, the model consists of a single 
two-state system coupled linearly to an infinite bath of harmonic oscillators, and described by 
the Hamiltonian 


H 


dk bibk do 


J(Ti — ha^ 


dk {bl + bk) <J3 , 


(1.4) 


where the Pauli matrixes and cxo = 1 describe the two-state system (“quantum spin”), b\ and 
bk are phonon creation and annihilation operators such that [bk,b\,] = k6{k — k'), [6^,6^'] = 
[^L ^I'] ~ t)are tunneling amplitude between the eigenstates of cxi is given by J, and h is 

an additional bias. The Ohmic dissipative two-state system is related to several other models, 
including the anisotropic Hondo model Ml], the resonant level model [12] and the inverse 
square Ising model [T2] . 

One important issue in the spin-boson model is the phonon-induced delocalized-localized 
transition. Such a delocalized transition at zero temperature is now considered as some kind 
of quantum phase transition. For the Hamiltonian fll.4l) the quantum transition of Kosterlitz- 
Thouless type occurs at = 1. The delocalized region 0 < < 1 corresponds to the antiferro¬ 

magnetic Hondo model, while the localized region corresponds to the ferromagnetic case. Here 
we will consider only the case 0 < 5 ^ < 1 . 

As it was pointed out in ref. [H] the Ohmic bath of harmonic oscillators can be interpreted 
as a simple bulk CFT - the massless Gaussian model. This allows one to reformulate the spin- 
boson model in the delocalized regime as a boundary flow problem, where the parameters of 
the Hamiltonian J and h play the role of the dimensionful boundary couplings. The flow starts 
from the Gaussian GFT with the Neumann (free) BG and with the decoupled spin degrees of 
freedom. For h = 0, the Gaussian held still satishes the Neumann BG in the IR hxed point; 
however, the spin proves to be completely screened (for details see, e.g., ref. | 1 ]). 

In this work we will study the vacuum overlaps fll.2p where the vacuums corresponds to 
different sets of the couplings (Ji, hi) and (J 2 , ^ 2 ). The arguments similar to that for the X-ray 
edge problem [IMZ] leads to the simple formula for the IR singularity exponent 


d2l = I ("12 


iriiY 


(1.5) 


where = (Gj |(J 3 | Gj)/ (Gj | Gj). Despite the lack of a rigorous proof, there are strong indica¬ 
tions, including numerical results from ref. [TS], that this is an exact relation for the spin-boson 
model with 0 < ^ < 1 . The aim of this work is to make steps towards the exact calculation of 
hdelities. 

The paper is organized as follows. In Sec. [2] we give a brief account of the basic concepts 
and facts and set up notations that will be used in the main body of the text. Sec. [3] reviews 
several well-known techniques for study of the orthogonality exponent and fidelities in the 
spin-boson model. The purpose of the next two sections is to develop a non-perturbative 
approach for a calculation of the fidelities. In the absence of IR divergences the Gell-Mann 
and Low theorem [19] allows one to express the vacuum overlaps in terms of the half-infinite 
time evolution operators in the interaction picture. However, the procedure which is based on 
the adiabatic switch of interaction generally fails for a system with gapless excitations. Our 
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approach is based on an axiomatic determination of the hdelities similar in philosophy to the 
form-factor bootstrap [20]. In Sec. 0] we argue that, in the case of spin-boson model, matrix 
components of the half-infinite time evolution operators can be interpreted as the quantum 
dost operators - the quantum counterpart of the dost functions for the pair of Sturm-Liouville 
equations. With this observation, the calculation of the hdelities can be considered as a part of 
the problem of quantizing the mKdV/sine-Gordon integrable hierarchy. The keystone element 
of quantum integrability is the Yang-Baxter type algebras with commutation relations dehned 
by certain quantum i?-matrix. In SecJS] basing on the results of the works [2T1|22], we propose 
a set of algebraic relations for the quantum dost operators, which is then translated into a 
set of functional equations imposed on the hdelities. Currently, the solution of the system of 
functional equations is known for the case hi = h .2 = 0 only. It was reported in ref. [23]. The 
last section of the paper contains a few remarks concerning the hdelities corresponding 

to a family of bare operators 

^ T ~ 

where (Xs G { 1, as, | (ai ± ia 2 ) } and a is a real parameter. Here we also present formulas for 
^i 2 (crl“^) in the case hi = h 2 = 0, which generalize the result of [23]. A derivation of these 
formulas is somewhat technical and it remained beyond the scope of this work. 


2 Preliminaries 

2.1 Basics of Gaussian model with Neumann BC 

We hrst consider the Gaussian model on the half-line whose dynamics is governed by the 
Hamiltonian, 

Iffree = ^ ^ dx + (4<I>)2 ) , (2.1) 

J—oo 

the Neumann BC, dx^{x^t)\x=o = 0, and the canonical commutation relations [*h(x), n(a:') ] = 
27ri g 6{x — x'), etc. The space of states splits up into the Fock spaces J^p - irreps of the algebra 
of creation-annihilation operators 

f2 dr 

J 4^ ^ = + ’ (2.2) 

whose highest weight vectors are defined by the conditions bk\p) = 0 {k > 0) and &o I P ) = 
^/^p\p). Throughout this paper, we will refer to \p) as p-vacuums. The Fock spaces are 
naturally equipped with the inner product dehned by the conjugation h\ = b-t and {p'\p) = 
3pf . 

Since Hfree = /o°° dk b^kbk, the Hamiltonian acts invariantly on each Fock space J^p. Further¬ 
more, all the p-vacuums correspond to the same zero-point energy, so that the ground state of 
the Gaussian theory with Neumann BC is a certain linear combination of the p-vacuums. With 
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Figure 2: The Wick rotation. In the Euclidean picture, the model can be interpreted as a 
two-dimensional classical statistical system on the half-plane ^e{z) < 0 with z = x + iy. 

the aim to dehne the ground state unambiguously, it is useful to consider the problem in the 
Euclidean picture where t is replaced by the Euclidean time via the Wick rotation t ^ y = t 
(see Fig.[2]). In the Euclidean picture all the helds are treated as functions of {x,y). Then, the 
ground state can be dehned through the asymptotic condition 

lim |vac) = |vac) , (2.3) 

y^—oo 

which holds true for an arbitrary real parameter a. The exponential fields act between the 
Fock spaces, : Tp i-4 J^p+a, and we will always assume the normalization condition 

{p + a \ Ip) = 1. Thus the ground state can be written in the form of a direct integral 

/ CX> 

dp I p) . (2.4) 

-CO 

The Hilbert space of the model is given by a linear span 

V. = span{ b-ki ■ ■ ■ b-k^ \ vac) | /cj > 0 & iV = 0, 1,... } . (2.5) 

Note that exponentials act invariantly on "H. 

The Gaussian model is manifestly invariant under the transformation ^[x,y) i—)■ —$(a:,p), 
which will be referred below to as C'-conjugation. The corresponding symmetry operator acts 
as 


C : Cbk = —bkC, CI vac) = I vac) . (2.6) 

Another global symmetry is the T-invariance. The antmnzfarp T-transformation acts according 
to the rule $(x,p) <F(x, —p) and the corresponding symmetry operator T can be dehned by 
the relations 


T ; T T = —, T | vac) = | vac) . 


(2.7) 


Finally let us note that the Gaussian held $(a:,p) splits into holomorphic and antiholomor- 
phic components 


<h(a:, p) = (j){x + ip) -f (j){—x + ip) . 


( 2 . 8 ) 
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In fact, the Gaussian CFT with Neumann BC can be interpreted as a model of a chiral bose 
field on the “unfolded” half-infinite line, whose Euclidean time evolution, (j){x,y) = (j){x + iy), 
is produced by the Hamiltonian 

Hiree = - - / dx (2.9) 

2vr^ J-oo 

through the commutation relation 

[(j){x2), ^(xi)] = ^ TTfif sgn(a;2 - Xi) . (2.10) 


2.2 Renormalization in the spin-boson model 

We now turn to the model of boundary interaction with the Hamiltonian 


H ^ CTo - (| ITg + /z) (Tg - J (Ti , 


( 2 . 11 ) 


where h/free has been defined by eq. fl2.1l) and H^ = n(a:)| 3 ,=o- This Hamiltonian acts in the 
tensor product of the Hilbert space fl2.5l) and the two-dimensional linear space whose endomor- 

phisms spanned by the conventional 2x2 Pauli matrices and do = f J ^ V The Hamiltonian is 


hermitian for real values of the parameters J and h. Without loss of generality one can assume 
that J > 0. In terms of the creation-annihilation operators fl2.2p the Hamiltonian H 

is given by eq. fll.4p . In this form it occurs as a particular realization of the Caldeira-Leggett 
Hamiltonian [5]. The model is usually referred to as the spin-boson model or dissipative two- 
level system and used to mimic the behavior of dissipative particle confined in a double-well 
potential. 

The spin-boson model needs renormalization. The Hamiltonian fll.4p has to be equipped 
with the ultraviolet cut-off A and consistent removal of the UV divergences requires the bare 
coupling constants J and g be given a dependence of the cut-off according to Renormalization 
Group (RG) flow equations. There exists a RG scheme where 


A dJ 


dA 


( 2 . 12 ) 


and because of this one can substitute the bare coupling J by the RG invariant energy scale 


E* 


const A 



(2.13) 


The latter is defined up to a multiplicative ^f-dependent constant and usually referred to as 
Hondo temperature in the context of the anisotropic Hondo model. The parameter h is in¬ 
terpreted as an external magnetic field applied to the impurity spin. It is often convenient to 
specify the Hondo temperature as 


E* 


dh^ 


Eo{J,h) 



h=0 


(2,14) 
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where Eq stands for the ground state energy considered as a function of the bare coupling J 
and h. 

Perhaps the simplest way to understand the renormalization scheme fl2.12p is based on an 
alternative form of the Hamiltonian H (12.lip . As it was already mentioned, the Gaussian 
theory with Neumann BC can be interpreted as a model of a free chiral boson. Eg. (12. 8 p 
implies that H^ = —2dx4>{0) and, therefore, the canonical transformation H HU with 

U = exp (icr30(O)) brings the Hamiltonian fl2.1ip to the form 

1 ^+oo 

H^ = U^HU = - / dx{dx(j)Y ao-has-+ (2.15) 

J —oo 


with cr± = |((Ti±i(j 2 ) and [cj){x), (j){x')] = ^ng sgia{x — x'). Notice that in eq. fl2.15p it is assumed 
that the Hamiltonian is expressed in terms of the renormalized exponential operators so that 
the bare coupling J = J{A) is substituted by the RG-invariant /i. In order to assign a precise 
meaning to the renormalized coupling one needs to specify a normalization condition for the 
renormalized exponentials. In fact, we have already accept the condition {p+a \ \p) = 1 . 

This sets a value of the leading term of the Euclidean operator product expansion 

g 2 i 0 (y) g- 2 i 0 (-y) ^ ^ {2y)-^3 as 1 / = it ^ 0 (l/ > 0 ) , (2.16) 

where Thus the renormalized coupling p has the dimension of [energy , 

i.e., /i = JA^ and p oc which yields formula fl2.13p . The exact E* — p relation can be 

extracted from the results of the Bethe ansatz solution of the anisotropic Kondo model [2^1125] 
(see also [ 22 ]): 

juvii + A 1 

g‘ = (i-g) ^ i f (r(i-j)^)---. (2.17) 

r(2 + 2(i?^) 


2.3 Fidelities A2\{0±) and ^421 (do, 3 ) 

Let us slightly generalize the setting from the introduction and consider the partition function 
Z 21 {p\y') of the half-infinite system with the insertion of a pair of hermitian conjugate boundary 
fields O and at the ends of the boundary segment where BG depends upon two energy scales 
(E^, ^ 2 ) (see Fig.[3|). In what follows we will make use of the notation 

Mo\v) = ^ ( 2 . 18 ) 

^11 

where Zn is the partition function of the system without any boundary insertions and whose 
BG is defined by {El, hi) homogeneously along the whole boundary. Of course, Z2i{0\y) is a 
complicated function and it is a challenging problem to compute it in a compact and manageable 
form even for the simplest boundary fields. However, since the physics of the model fl2.1ip is 
well understood now, one can make some general predictions regarding its behavior for small 
and large values of y. In what follows we will mostly discuss the simplest case with O given by 
the diagonal matrixes: cxo, (J 3 , or 

0± = ^ (aoEas) . (2.19) 
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Figure 3: An universal part of |(G 2 | C^(0) I )| can be extracted from the partition function 
of the system with inhomogeneous boundary, and with the insertion of a pair of hermitian 
conjugate boundary helds O and Ol 


Since the model is asymptotically free at short distances, the effect of the boundary energy 
scales (i ?2 5 £- 2 ) from the interval (0, y) becomes negligible wherein its size shrinks to zero, and 
therefore 


Here we use 


\i-m Z2i{0\y) 


(1 ± mi) 


for O = do, 0-3 
for 0 = 0± 


m = 


{G I (T 3 1 G) 
(G|G) 


( 2 . 20 ) 


( 2 , 21 ) 


and its subscript “ 1 ” shows that the expectation value is taken over the ground state | Gi) 
corresponding to {El hi). Note that m coincides with the impurity magnetization in the 
context of the anisotropic Kondo model. It can be expressed in terms of the exact ground state 
energy Eq of the Hamiltonian fl2.11l) as 


m = — 



( 2 . 22 ) 


The spin degrees of freedom are freezing at large distances so that the large-?/ asymptotic of 
Z 2 i{ 0 \y) has a form similar to fll.ip with the same exponent ^21 = <^21 ({£'*, hi}) for any choice 
of the diagonal matrix O E (cxo, < 73 , 0±}: 

Z 2 i(Oly) = \A2 i{0)\^ (l + o(l)) as ^ +00 , (2.23) 

where AE 21 = Eo{Elh 2 ) — Eo{Elhi). Although (I2.23p dehnes A2i{0±) and A 2 i(cro, 3 ) in 
absolute value, their relative phases are dictated by the relations 


A2i(<7s) — ^21(0+) -f (—1)^ A2 i((9_) (s — 0,3) . (2.24) 

In what follows, we will call A 21 {O) as “hdelities” and treat them as scaling functions depending 
on the magnetic moments mi, m 2 fl2.22p as well as the pair of the Kondo temperatures El 
E 2 (12.141) . In fact, since they are dimensionful quantities, they depend non-trivially upon the 
dimensionless variables mi, m 2 and a = log(F' 2 /-£'i) only. 





















3 Mean field, perturbation theory and Toulouse limit 

In this section we ontline several common approaches for study of the hdelities in the spin-boson 
model. 

3.1 Mean field approximation 

We start with the model which is considerably simpler than the spin-boson model. The sim- 
plihed Hamiltonian is obtained from H (I2.1ip through the substitution of the Pauli matrix 
by a constant m. It splits into two non-interacting parts: 

TTi 

= (3.1) 

and a 2 X 2 matrix H2 = —h cxs — Jcti. This can be thought of as a mean held approximation, 
with the value of m = m{h) given by the relation m = — where Eq = —y/J'^ + h'^ is the 
lowest eigenvalue of H2- 

To construct the vacuum state for the mean held Hamiltonian one can use the interaction 
picture with the term oc m in fl3.ip is treated as an interaction. The unitary operator U{t, to) = 
giHfree(i-io) calculated explicitly in this case: 

U{t,to) = ^ 

where <I)(0, 0) 'ppg vacuum for (13.ip is obtained through the Euclidean 

time evolution of the state |vac) fl2.4p . Thus the vacuum state for the whole mean held 
Hamiltonian is given by 

SJ > = eT *»<“' hm e-T *»“«> | vac) ® . (3.3) 

to^oo ' ' y/2 Vvl -m) ' 

where the limit is taken along imaginary time direction yo = ito ^ —oo as it is shown in Fig.[2l 
Taking into account the dehning property fl2.3p of the state | vac), one obtains 

I ) = r+ I vac) (g) I t) -F T_ I vac) (g) | |) , (3.4) 


where 


T+ = 


1 ± m 


e 2 


if- i>b(o) 


(3.5) 


and we use the common notation for as-eigenvectors. Similarly one has 

{ hi I = (vac I Q+ (g) (t I -f- (vac | Q_ (g) (I | 

with 


= (r±)'= e-*?**™ 


(3.6) 
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It is easy to see now that the mean held approximation yields the relation 

Z 2 i{ 0 ±\y) oc (vac I | ^ = 1 ( 3 . 7 ) 

and therefore to eq. fll.5p for the IR singularity exponent. As for the hdelities A2i(0±'), it is 
worth to note the relation 

(p'|gi'^Ti'^|p) = A2l((9±) 62 p'- 2 p,m,-m 2 • (3-8) 

Within the mean held approximation A2i(0±) = ± m 2 )(l ± mi), which is found to be 

an adequate approximation as —?■ 0. Contrary to the hdelities, the formula <^21 = f 
is expected to be exact as mi and m2 are understood as vacuum expectation values of as in 
the spin-boson model. Evidences in its favor are presented in the next two subsections. 

3.2 Renormalized perturbation theory 

Similarly to the case hi = h 2 = 0 considered in ref. [23], the hdelities A2i{0±) with hi, ^2 7 ^ 0 
can be calculated by means of the renormalized perturbation theory in coupling g for the 
Hamiltonian fll.dp . The result of perturbative calculations turn out to be in agreement with 
the orthogonality exponent ^21 = \{fn 2 — miY. This ensure that A2i(0±) can be written in 
the form 

A,_j(Og = _ (3.9) 

where 

a = log , (3.10) 

and the prefactor has a dimension of [energy, so that A± are dimensionless amplitudes. 
Notice that in the case /i 2 = hi, ^2 = hi, the orthogonality exponent vanishes and A± should 
satisfy an exact relation 

hl±(m, m I 0 ) = ^ (1 ± m) . (3-11) 

It can be shown with somewhat cumbersome but straightforward effort that 

hl+(m2,mi I «) = ^ V(1 + rn2){l + mi) (2e'^-®) (3-12) 

X (1 - ml) f d-”^ 2 )(mi-m 2 ) I(l-mi)(m 2 -mi) 

X ^1 + I ^ (1 - miY + (1 - m 2 )^ - (1 - mi){l - m2) S coth ^ + 0{g‘^) ^ 

and 

hi-(m 2 , mi I a) = hl+(—m 2 , —mi | a) . (3.13) 

Here '^e stands for the Euler constant, and we use 5 = a + 7 log The quoted result 

shows that the perturbative amplitudes A± are multivalued functions of the complex variables 
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(m 2 , mi I a). However their overall phases can be chosen in such a way that they are real 
analytic within the domain 

D = I (mi, m 2 I a) : mi ^2 G (—1,1) ; a G M } . (3-14) 

It is also easy to see that A± satisfy the condition 

A±{m2,mi\a) = A±{mi,m2\ — a) . (3.15) 

Besides, the use of the perturbation theory allows one to determine the relation between m 
and dimensionless ratio h/E* in a form of a power series in g. In particular, to the first order, 
one has 

I? = 7t= 7P (^■ f■ P-'®' 

With this formula the perturbative amplitudes can be expressed in terms of hi/E* {i = 1,2). 
Notice that, from the Bethe ansatz solution of the anisotropic Kondo model, it is known that 

[2ll[25] 


i 



den A 1-9 


M{u) 
cn + iO ’ 


M{u) 


r(i - + t) 

x/?r(i-^) 


(r(i-9)) 


(3.17) 


where 


1 


1-9 

r(i-9) 

[v?(i-9)r(i + ^)J 

\ h ) 


(3,18) 


This remarkable exact result implies the following general structure of the perturbative expan¬ 
sion fl3.16p : 


j / 00 n—1 \ 

— = ( 1 + m^M , (3.19) 

n=2 1=1 ' 

where Cni are some numerical coefficients. 

3.3 Toulouse limit 

In the case e? = |, which is sometimes referred as to “Toulouse limit”, the Hamiltonian 
fl2.15l) can be fermionized in terms of the chiral complex fermion field 

^^Toui = —: / dx dxij — g-[d)ij{^)-\-+ h {(Ad — dd)) ^ (3.20) 

27ri J_oo 

where {'0^(a:),'0(a;')} = 5{x — x'), {d^,d} = 1, etc. Equivalently, the model can be understood 
as a boundary flow for the Dirac fermion, massless in the bulk. In order to construct an 
Euclidean action governing this boundary flow, we define 'ijj{x,y) = 'i/j{—x,y), so that the 
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Hamiltonian fl3.20p for // = /i = 0 takes the form ^ dx — 'ip’^dx'ip) ■ The helds 

^|J and '0 are interpreted now as components of the massless Dirac fermion, both dehned on 
the half-inhnite lane x < 0 and satisfying the bulk equations of motion = dz'ip = 0 with 
z = x+iy, z = x—iy. The complex fermions can be substituted by two Majorana-Weyl fermions: 
^ ^ ^ ^ Qf j,ggi fermions {ipj, ipj) satishes the free BC, |a;=o = 0 

and, as it was explained in ref. [26], should be described by means of the action %] 

which involves additional boundary fermionic degree of freedom aj = aj{y): 




dy 


- - 1 f 1 - \ 

dx +- dy{^—{'ijjj'iljj)\x=o + ajdyajj. (3.21) 


An Euclidean action, corresponding to the Hamiltonian fl3.20p with non-vanishing couplings y 
and h, is given by IZS0 


where 


2 

•TtouI = •TMwlV'j) V'j) + 
i=i 



(^y&Biy) - h Msiy) ) , 


2 

©B = - ^ ttj [ijj -f ipj) |3,=o , = 2i OiOa . 

i=i 


(3.22) 


(3.23) 


Clearly, the Grassmannian boundary helds d = d* = in the path integral formalism 

correspond to the nilpotent operators d and d^ in the Hamiltonian picture. 

The general solution of the bulk equations, dz^J = dzip = 0, are given by the Fourier integrals 

/ CO poo 

dk Ck e‘^^ , V’ = / dk Ck . (3.24) 

-OO j —CO 


In their turn, the boundary equations of motion corresponding to the action 03.221) allows one 
to express the Fourier modes Ck and the Heisenberg operator d{y) in terms of c^: 

= ^ith 

k + my? -2h ^ ^ 

and 


d{y) = i y 



Ck e 

k -|- iTTyU^ — 2h 


(3.26) 


With this, it is straightforward to compute the Euclidean propagator for the complex boundary 
fermions: 


( d*{y) d{0)) = D{y\- h) Q{y) - Q{-y) D{y\h) , (3.27) 

^ Although in ref. |2^ the case of a single Majorana-Weyl fermion was only discussed, Eg. (13.221) is an apparent 
consequence of the Chatterjee-Zamolodchikov result. 
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where Q{y) = |(1 + sign(|/)) and 


l+m 


D{y\h) = /i^ / dk 


^-My\ 


{k — 2hy + TT^/i^ 


-)■ 


l+m 

cos^(^) 

AE*y 


In the last formula we use the notations E* = 


as 2/ —t 0 
as y ^ oo 


(3.28) 


2 f irh 

m = — arctan - 

TT \2E* 


(3.29) 


which are consistent with general relations fl2.17p and fl3.17p taken aX g = \. 

In order to calculate the orthogonality exponent, one needs to introduce an explicit IR 
regularization. Let us restrict values of the space coordinate x to the segment [—L, 0] and choose 
the free BC at x = —L\ (-0 — 'ip)\x=-L = or equivalently, Ck = Ck e“‘^^. 

Taking this together with eq. (l3.25p . one obtains the quantization condition 2Lkn + 2 6{kn) = 
27m {n G Z), so that Fourier integral expansions (I3.24p should be substituted by discrete 
sums '0 = 'YYY=-oo^ri and similar for Y- Note that the real-time evolution of the 

fermion modes c„ is produced by the Hamiltonian H = Yln=-oo through the canonical 

commutation relations {cl,c„/} = Sn^n', = (c^)^ = 0. A ground state of the system of 
fermionic oscillators is dehned by the requirement that all energy levels bellow the Fermi level 
kp = 0 are occupied. In this situation, according to Anderson [2], the orthogonality exponent 
for the ground state overlap (112 I ) is determined by the the difference of the phase shifts at 
the Fermi level; 


<^21 



^ (<5.(0)-i,(0))V 


(3.30) 


In its turn, the phase shift 6{k) fl3.25p at fc = 0 can be written in terms of m fl3.29p as 
5(0) = I (1 — m), and therefore ^21 = | ("^2 — miY. This coincides with fll.Sp specialized at 
9 = i 

It is useful to note that the orthogonality exponent can be written in the form 


^21 — 2 (^'i ~ 92 Y , 


(3.31) 


where qi = (Hj | g | ) are vacuum expectation values of the operator q = Nu¬ 

merical results presented in ref. |IH] suggest that eq. (l3.3ip remains valid for g 7^ For their 
calculations, the authors used the resonant level model na whose Hamiltonian is obtained by 
adding a four-fermion interaction term to Hpoui (I3.20p . 

F/rlm = HtouI + U {d^d - (Yd) : V'V(O) : • (3.32) 


As it is well known (see, e.g., [271128]), this is a fermionic version of the Hamiltonian for 
general values of g G [0,1]. The difference qi — q 2 in (13.31 p was referred in ref. [H] as to the 
“displaced charge”, which is, in a sense, the difference between vacuum expectation values of 
the charges associated with the global U{1) symmetry of the resonant level model. We may 
now note that formulae fll.Sp and fl3.3ip for the orthogonality exponent coincides provided 


qi-q2 = 



{m2 - mi) . 


(3.33) 


13 





























The last relation indeed holds true for the resonant level model (see Chapter 28IV.2 in ref. ^7] 
for details). 

Returning to the spin-boson model at the Toulouse limit, we acknowledge the relations 
between the normalized partition functions = Z 2 i{az,i\y)\E*=Ei,h 2 =hi two-point 

functions of the boundary helds (I3.23p : 


Zii{(y^\y)\, = {mB{y)mB{Q)) (3.34) 

'y —2 

Zii{(Ti\y)\g^i_ = (6b(i/)© i?(o)) . 


The boundary equations of motion corresponding to the action (I3.22p . allows one to represent 
the Heisenberg operators 9 JIb(|/) and ©s(|/) in the form of normal-ordered combinations of d{y) 
and its Hermitian conjugates: 


TIb = 2 : dd^ : +ni 

&B = (^{dyd) d^ + {dyd^)d — AhSd) +{(Sb) , 


where the vacuum expectation value, 


( 6 b) 


TT 


2 log(f) +log(l+ 




(3.35) 


(3.36) 


contains the non-universal term which depends on the UV cutoff scale A. Using the Wick 
theorem, the two-point functions fl3.34l) can be written in terms of D± = D{y\ ± h) fl3.28p . 

{mB{y)MB{0)) = m‘^+ 4D^D+ (3.37) 


D~D+ (D~ 7)+ 

( GB{y) Gb( 0) )conn = 


D_ 




Here the dot stands for the derivative w.r.t. the Euclidean time y and the abbreviation “conn” 
means the connected correlation function: ( ©^(y) ©b(0 ) )conn = ( ©s( 2 /) ©b(0 ) ) — (©s)^. 
Similarly, one has 


/i \D+ 



(3.38) 


Since the two-point functions (I3.37p . (I3.38p are available in closed forms, the partition 
function 22 i(o'o|l/) can be calculated perturbatively in powers of 6 m = m 2 — mi 1 and 
6 a = {E 2 — E\)/E^ 1. In the case 6 m = m = 0, details of the calculations may be found in 

Appendix in ref. [23]. Similar calculations for non-zero m and 6 m yield the first non-vanishing 
terms of the Taylor expansion of the fidelity A 2 i(cro) in the Toulouse limit: 


Mi{cro)\ 1 — 


TTCOS 


Ae<EE* 


d2i 


1 - Ac 


((5a) 


- A»i 


{ 6 m) 


+ 0{6^ 


(3.39) 


Here m = mi, E* = Ei, 621 = {6m)‘^/16 and 

Aoa 


1 / TT , , 

—^ 1 H—(m -I- 1) tan , , 

47r2 V 2^ ^ V 2 7 


Ac 


TT" Aaa 


sm 7rm 


cosi^vrmj 

16sin2(^) 


^7rm j 
1 -|- Tim tan 


TT , , / TTm \ 

l + -(m-l) tan(^—j 


(f) 


(3.40) 
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4 Jost operators 

In this section we introduce the notion of quantum Jost operators and discuss their properties. 


4.1 Classical Jost functions 

Let us hrst consider the limit g ^ Q where the boundary held ns(t) = is treated as 

a classical held. In this approximation the Hamiltonian fl2.1ip describes the quantum spin in 
the presence of a classical background held. The corresponding non-stationary Schrodinger 
equation has the form 

=-(Jo- 1+ 0c(t) 0-3) T , (4.1) 

where we use | -|- ht, satisfying the asymptotic condition 

4 >c{t) ht + 0(1) as t ^ —00 (4.2) 

with h > 0. Let U{t,to) be a time evolution matrix for fl4.ip . It can be written as 

U{t, to) = S{t, to) , (4.3) 


where S stands for a time-ordered matrix exponential 

S{t, to) = 71 exp (^i ^ dt J (ct_ + ) ) • (4-4) 

Consider the time evolution of the spin-up state | f) starting from the initial moment to- The 
following limiting behavior for to — 00 , can be easily established: 


t) ^ t) + T e-'C-'')'” (4.5) 

+ ( ^«'-!(«) 11) - 4 u >) . 

Here '^±^{t), T±]5(f) are the Jost solutions of the Sturm-Liouville equations 

{-d^ + U{t)) , U = + i 0, (4.6) 

{-d^ + U{t)) T±, = U = 

satisfying the asymptotic conditions at t ^ — 00: 

4'±k(t) ^ , T±,(t) ^ e±'"* , (4.7) 

where k = \/+ h‘^ > 0. Because of the presence of oscillating phase factors, the r.h.s. of 
fl4.5p does not possess a hnite limit as to tends to —00. However, if we assume that the coupling 
J is switched oh adiabatically. 


lim J{to) = 0 , (4.8) 

tQ—^—00 
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then limi(,^_oo k = h > 0 and the hrst term in fl4.5p has a hnite limit. The second term will still 
oscillate oc as to — cxo. With the aim to suppress these oscillations, let us fix the value 
of t, say t = 0, and assume that the asymptotic behavior fl4.5p holds true in the infinitesimal 
wedge 0 < arg(—to) < e —)> +0 of complex plane to- Then, taking the limit along any ray inside 
the wedge, one obtains 


arg(-tQ) = + 0 V 



where we use m 


■ The t-independent connection coefficients 


(4.9) 


rf = '■'+..(0). rt' = 4,^,(0) (4.10) 

are sometimes referred as to dost functions. It is easy to see that eqs. fl4.6p . fl4.7p imply that 

T±k(t) = {dt - i0c)T±k(t) , T±k(t) = - {9t + i0c)T±k(t) 

and also = \h^]j(t) for real J, h and t. Using this and also taking into account that the 

Wronskian of T_]5.(t) and equals 2ik, one finds the bilinear relation 

Comparing fl4.10p with eqs. fl3.5p . fl3.6p we note that within the mean field approximation, the 
classical connection coefficients T_|_]j(0) and are substituted by the exponential operators 

e^'^B(O) and e“^'*’®(0), respectively. 


4.2 Anticipated properties of Jost operators 

Motivated by the above consideration, we start from the Hamiltonian fl2.15p which describes 
an interaction of a local spin impurity with a chiral bose field on the whole line. Consider the 
interaction picture, treating the term oc /i as a perturbation. Let us perform the Wick rotation 
from the very beginning, so that (j){x, y) = (j){x + iy) for the chiral bose field in the interaction 
picture and introduce the ^/-ordered exponent 

S{y 2 , yi I /iL, h) = Ty exp dy miy) e'^^^a- + (T+) ^ . (4.12) 

Here the renormalized coupling p is substituted by |/-dependent function, fiiiu) > 0 such that 
h = /^l( 0)) which is switching off adiabatically within an Euclidean time interval \y\ < L. 
For technical reason, it is also convenient to choose the smoothing function to be an even 
function of y. The explicit IR regularization allows one to define 

U^-\pL,h) = 5(0,-cx) I/iL, h) (4.13) 

U^+\fiL,h) = 5(+oo,0|/Xi,/^)e-‘‘^(°’°)"^ 
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which are operators acting on the impurity spin, whose matrix elements are themselves oper¬ 
ators acting on the free bosonic degrees of freedom, i.e. in the Hilbert space "H fl2.5p . In the 
absence of interaction and for h > 0, the vacuum state is given by the product | vac) 0 | t )• 
Naively, the vacuum in the interacting theory occurs through the large-L limit: 


eA/(L) 


vac) 0 1 t) 

— y 

\VL) 

(4.14) 



L —^~l“Oo 

eA/R) 

(vac 0 (t 


— y 

L —^-|-(30 

(H| , 



where A/(L) = J^^dy [Eo(^^L(y), h) -|- h). But, because the IR problem, these asymptotic 
relations cannot be literally true. We will try to overcome this obstacle using a heuristic picture 
which is based on the notion of quantum dost operators. Namely, we postulate that the exact 
vacuum state is given by relations similar to the one obtained in Sec. l3. II within the mean field 
approximation, i.e., 

|H) = T+(/i, h) I vac) 0 I t)-f T_(/i, h) I vac) 0 I I) (4.15) 

(H| = (vac I Q+(/i, h) 0 (t I-f (vac I Q_(/i, h) 0 (I I , 
and Tg and Q^, with e = ± and y,h > 0, are operators acting as 

h) : H ^ H Sz J-p 

Qei^y, h) : "H "H & J-p e-^ i (4.16) 

such that p-vacuum expectation values involving bilinear combinations of T and Q are expressed 
in terms of the dimensionless amplitudes A± (13.9^ : 

{p' \Q±{^J'2,h2)T±{yl,hl) Ip) = A±{m2, mi | a) 5 

2p^—2p,mi—m2 * (4.17) 

Notice that, since T and Q act invariantly in JL, the p-dependence appears here through the 
Kronecker delta only. In what follows T- and Q-operators are referred as to dost operators. In 
order to predict their properties, we shall invoke the intuition which is based on the classical 
limit, the results of perturbative calculations and global symmetries of the model. 

• T — invariance. The time reversal transformation acts as e-)■ and 

a± (Tip, (Ts e-)■ —(T3. The antiunitary operator defined by fl2.7p . satisfies the relation 

gia20(O,y2) . . . gian0i(O,j/„)jrjp _ gia„0(O,-j/„) _ _ _ gia20(O,-j/2) g-iai0i(O,-yi) 

for Pi > P2 > • • • > Pn- Ignoring the problem with the IR divergency, one can expand the 
p-ordered exponent (14.121) and find 

TTg(/i,h)T = gg(/i*,h) , Tgg(/i,h)T = Tg(/i*,h) . (4.18) 

Notice that the T-invariance, when it is applied to (I4.17p . leads to the relation (I3.15p . 

• H ermiticity . The dost operators satisfy a hermitian conjugation condition 

Qe{p, h) = (Tg(/i, h))^ (/i, h > 0) , (4.19) 

which follows from the relations (T_)^ = (T+. 

Toffether with fl3.15p. the coniuffation implies the realitv condition iAA'mo.m,^ l«))* = 
A±{m 2 ,rni\a) for real positive /i*, hi {i = 1,2). 
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• C — invariance. The formal C'-transformation acts as g-2ia</)(o,y)^ 

(Js 1-^ —(Ts- It is well known that there is no spontaneous magnetization in the anisotropic 
Kondo model. Because of this we expect that there exists an operator C such that = 1, 
and the dost operators for h < 0 can be introduced through the relations 

T_,(/i,-h) = CT,(/i,h)C (4.20) 

Q-e^,-h) = CQe{fi,h)C. 

An important consequence of the Hermiticity, C- and T-invariance is that p-vacuum 
expectation values of dost operators are expressed in terms of a single, real analytic 
function of m G (—1, 1)0 

{p'\T±{fi,h)\p) = F{±m) 62 p'- 2 p,m (4.21) 

{p' \Q±{p,h)\p) = F{±m) S2p-2p'-m ■ 

As for the matrix elements fl4.17p . the C-invariance implies that A± satisfy the relation 
fl3.13p within the principal domain fl3.14p . 

• Bilinear relation. As it follows from eq. (l3.1ip . ^+(m2,mi |0) +^_(m2,mi |0) = 1. 
Taking this into account along with the conjugation rule (I4.19p . the quantized version of 
the bilinear relation (14.111) is expected to be given by 

J2Qe{Fh)T,{p,h) = l\n . (4.22) 

€=:L 

• Lorentz invarianee. Let us introduce the complex coordinate 2; = x + iy and the 
polar angle ^|J = a.Tg{iz). The generator of infinitesimal Euclidean rotations coincides 
with {—iK), where K is the Lorentz boost generator. It is crucial for our analysis that 
the angular evolution of dost operators turns out to be remarkably simple. Namely, 

(4.23) 

Qe{fi,h) p,e"^h) . 

This follows from three observations; First, the exponential helds in (14.121) are chiral (holo- 
morphic) = e^^“^(z)) with the Lorentz spin g. Second, the exponential operator 

located at the origin is not affected by the Euclidean rotation; g±i0(o,o) _ 

g±i(/>(o,o)_ have to accept that integrals containing combinations of the holomor- 

phic fields are not changed by rotations of the integration contour in the limit L —)■ cx). 

The simple geometry of the Euclidean plane suggests that T- and Q- operators are related 
through the Euclidean rotation of angle tt, combined with the C-conjugation; 


Q,(/i,h) = Ce-'"'^T,(p,/i) e'^'^C, 

(4.24) 

or. eauivalentlv (see eas.(l4.20ll.(l4.23l)j 


g,(p,h) = T_,(e^"(i-^V,/^) • 

(4.25) 


^ Notice that as it follows from (14.101) . limg_,.o F{m) = {^-^) ^ ■ 
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Figure 4: The contour of analytic continuation in eq. fl4.26p as mi G (—1,1) and a G M. 

It is worth to keep in mind that these formulae should be understood in a weak sense 
as relations for the analytic continuation of a certain class of matrix elements of Jost 
operators. The applicability of eqs. fl4.23p (which is crucial for deriving fl4.25p l requires 
that the contours of integration, involving in the construction of Jost operators, can be 
rotate freely within the Euclidean plane. This may be not the case for general matrix 
elements. 

We expect that fl4.25p can be applied for the matrix elements fl4.17p . and therefore the 
dimensionless amplitudes A± are related by the analytic continuation with two-point 
functions containing T-operators only: 

{p'\T^{H2,h2)T±{ni,hi)\p) =A±{e~^""m2, mi\a-m) J2p'-2p,mi+m2 • (4-26) 

Here the phase rotation e“'’^ means the analytic continuation along the clockwise half¬ 
circle of radius smaller then one (see Fig.HP for —1 < mi < 1 and a G M. Notice that the 
analytic continuation appearing in the r.h.s. of fl4.26p does not involve any problem for 
the perturbative amplitudes (13.121) . fl3.13p FI This supports rather sweeping assumptions 
that have been made in derivation of 04.261) . 


5 Algebra of Jost operators 

At the moment it is not clear how to deal with the hdelities beyond the scope of perturbation 
theory for arbitrary values /ij and hi. However, in the case oi hi = h 2 = h much can be said 
about the matrix elements 

{p \ T^.,{p2,h)Te^{pi,h) \ p) 

~ ^£2^1 (m2, mi) 5 

2p'—2p,m2-\-7ni • ( 5 . 1 ) 

Notice that according to eq. 04.261) 

^±(m2, mi = F^±{e'"^m2,mi) , (5.2) 

^ In this connection it deserves mentioning that ea. (l3.13p should be understood as Ae{m 2 ,'mi\a) = 
Aj(e^*'^m 2 , e^''^mi | a), which is satished for any sign ±. 
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and a is not an independent variable as hi = /i 2 - In this case a can be written as 


“ = (^) ■ (5,3) 

where A* [i = 1, 2) stands for the dimensionless ratio which is a certain function of m*, 

i.e., Aj = The corresponding inverse function m = m(A) is given by eq. fl3.17p . In the 

present discussion we will use both variables and Aj assuming that they are related through 
the formula fl3.17p . Since h > 0 is assumed to be hxed from now on, it makes sense to simplify 
the notation for the dost operators: 

h) = Te{m) , h) = Qeim) . (5.4) 


5.1 Commutation relations 

Let us recall that y) ai y = 0 satisfies the commutation relation fl2.10p . Since the Euclidean 
time dependence of the chiral field y) is very simple, we can translate fl2.10p into the 
commutation relation at x = 0: 

[0(0, 1 / 2 ), 0(0, 1 / 1 )] = ^ sgn (|/2 - yi) . (5.5) 

According to ref. [22], if the matrix-valued operators 2 /i I hi, h) {i = 1,2), 

where S is defined by eq. fl4.12p with 1/2 > yi and with /i* are set to be ^/-independent constants, 
then the Yang-Baxter equation is satisfied in the form 

R{Li®l){l®L2) = {l®L2){Li®l)R. (5.6) 


Nontrivial matrix elements of the 4x4 matrix R read explicitly 


Rtt = RZZ = l, 


RtZ = RZl = 


A? 


Ai 


qA) 


q-Ui ’ 


RZ+ = R-^ 


qAf-q ^A| 


where we use q = e'’^^ and Aj = Although a mathematically satisfactory construction 

of the dost operators is absent, the arguments similar to those from ref. [2T] suggest that the 
dost operators satisfy the commutation relations 




£l 




(5.8) 


£i,4=± 


where the i?-matrix obeys the Yang-Baxter equation together with the “unitarity” and “crossing 
symmetry” relations 


£'i,£2=± 


R<l'Rmi,m2)RTJ(m2,mi) = , 


(5.9) 


£'j,£'2,e^ = ± 


^£i+£',0 ^2) i?^g(mi, m2) = 
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The formal role of the Yang-Baxter and unitarity relations is clear; the unitarity is reqnired 
for self-consistency of flS.Sp . whereas the Yang-Baxter equation is the associativity constraint. 
In its turn, the crossing symmetry allows one to supplement algebraic relations flS.Sp with an 
extra bilinear relation, 

T^{m) = 1 , (5.10) 

e=:L 

which follows from eqs. fl4.22p and fl4.25p . 

The i?-matrix in the form 


= R{mi,m2) , (5.11) 

meet all the necessary requirements, if the normalization factor satisfies the conditions 

i?(mi, m 2 ) i?(m2, mi) = 1 (5.12) 

a\2 _ Q-i\2 

i?(mi,m2) i2(e'’"mi,m2) = —k- \T^ ■ 


(recall that we use the convention A* = A(mj)). Some extra conditions are imposed by the 
global symmetries. Namely, the Hermiticity and T-invariance require that 


F(mi, m2)F*(mi, m 2 ) = 1 

whereas the U-symmetry yields the relation 

for 0 < mi ^2 < 1 , 

(5,13) 

F(e^'"mi,e^'"m2) 

= F(mi,m2) . 

(5.14) 


5.2 Normalization factor 

The algebra of the dost operators together with the global symmetry relations and normaliza¬ 
tions conditions lead to a system of functional equations imposed on the two-point function 
F£ 2 £^(m 2 ,mi) fl5.ip . Assuming that 0 < mi ^2 < 1, the system reads as follows: 


e'j,£'2=± 


(5.15) 


and 


Fe2£i(e'’"m2,mi) = F.^^ _£2(e‘’'mi,m 2 ) = F_e2-£i(m2,e '’"mi) = (F£2£i(e‘’"m2,mi))* (5.16) 
and 

F_££(e‘’"m, m I 0 ) = ^ (1 -I- em) . (5-17) 

Using the perturbative results fl3.12p . one can check that all the conditions for F=p± are satisfied 
at the first perturbative order, provided the F-matrix has the form fl5.1ip with 

logR{mi,m2) = ^9 (1m2mi)-F 0(/) . (5.18) 

2 m 2 -|- mi 
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Without additional analytical assumptions the set fl5.12p - fl5.14p and fl5.18p does not unam¬ 
biguously dehne R{mi,m 2 )- However, the hrst perturbative correction allows one to make a 
conjecture about the exact normalization factor. Namely, a simple calculation shows that, at 
least at the hrst perturbative order, the normalization factor obeys the condition 


a d 

^2 Tv;- 1 “ -^1 


dX, 


dXi 




17111712 


(5.19) 


Here = m{Xi) {i = 1,2). If we accept fl5.19l) for any 0 < (7 < 1, the reconstruction 
of logi?(mi,m 2 ) requires a minimum amount of additional analytical assumptions. Indeed, 
consider the limit when h —)■ 0, keeping the Kondo temperatures and E 2 hxed. In this case 
both mi and m 2 turns to be zero, but their ratio remains hnite. Then 


hmi?([7e“b [7e“2) = Ro{ai - 02 ) , 

where ai — a 2 = log^E^/El), and as it follows from eqs. fl5.12p . 

sinh(l — g) {a + ivr) 


Ro{a)Ro{—a) = 1 , Ro{a + i7r)Ro{a) = 


sinh(l — g)a 


(5.20) 


(5.21) 


(Notice that for unrelated hi and ^2 the variable a coincides with the one dehned by eq. fl3.10p .) 
The “minimal” solution of (15.2 ip (i.e., such that logi?o(o) is analytic in the strip 0 < '^i7i{a) < 
TT and bounded at a —)■ -|-oo) has the form 


i?o(«) = exp —i 


do; sin(Q;a;) si'^h 2 (i-g) 
a; cosh^ 


(5.22) 


Eq. fl5.19p . as a linear partial differential equation subject of the asymptotic condition fl5.20p . 
can be easily solved by means of the Fourier transform. Using fl3.17p one hnds 


R{i7ii, m 2 ) = exp 
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00 


CJ2 - 


iO 


(5.23) 


(wi-h iO)(a;2iO) 002 + ooi 

The function fl5.23p meets the required conditions fl5.12p - fl5.14p . It also possesses a small -(7 
series expansion 


logi?(mi, m 2 ) = ivr E g" f„(mi,m 2 ) , 


(5.24) 


n=l 


with 


1 m 2 - mi 

n = --^- (1 + mi m 2 ) 

2 mi -|- m 2 

mim2(l-m?)(l-m^) f ^ 2,0 1 

72 = -- 2 a 2 - mi - m 2 2 mim 2 log — 

(mf — m^) V Vmi// 

mim2(l — mi)(l — m^) 


(5.25) 


rs = 

12 


m; 


m: 


2A3 


(mj — m^) (1 -|- mj -|- m^ — 3mim2) 


(mi — 1712^ {ml -|- m^ -|- 17111712 — 1) -|- 2 I7iii7i2{2i7il ml — ml — ml) log^ j ^ ^ etc. 
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It would be valuable to confirni these higher-order corrections using the renormalized perturba¬ 
tion theory. More importantly, the origin of fl5.19p (if it is true, of course) needs to be clarified. 
It may be useful to note here that the relation fl5.19p can be equivalently written as 

( d d \ ~ f d d \ 

V ^ ^ "^ 2 ) = ITT A 2 ^ - Ai — j d2i{-m2, mi) , (5.26) 

where .R(mi, m2) = i?(mi, m2) and ^21 = f (m2—mi)^ is the orthogonality exponent. 


5.3 Fidelities ^21 ((^0,3) for h = 0 

As h —)■ 0, the system of functional equations fl5.15p - fl5.17l) is simplified dramatically. In this 
case 


lim F+zr 


([)e“’,l)e”‘) = /(oi - 0 : 2 ) 


(5.27) 


and / satisfy the equations: 


f{i7r + a) = /(ivT-a) , 


/(-«) 

/(«) 


sinh (ivr -|- a) 

sinh (Itt — a) 


(5.28) 


Let us assume now that log / is analytic in the strip 0 < ^m{a) < vr and grows slower then 
any exponential of a at inhnity0 With these analytical assumptions, the functional relations 
fl5.28l) define / up to a constant multiplicative factor. Then the normalization condition fl5.17p 
hxes / unambiguously. Finally, using relations fl5.2p . we derive 


^ 21 ( 0 - 0 ) I/ii=h 2 =o = 2/(a - i7r) = «(«) 

^2l(o'3)Ui=/i2=0 = 0 , 


sinh(l - ^) f 
(1 - g) sinh f 


where 


u{a) = exp 


irgu) 

afWTy 


do; sin^(Q;Ci;/2) sinh 
oj sinh nu cosh ^ sinh 2 (^-g) 


A numerical verihcation of the prediction fl5.29p was performed in ref. 


(5.29) 


(5.30) 


6 Concluding remarks 

Although the problem of finding hdelities A 2 i{(Ts) for non-vanishing h remains unsolved, it 
seems useful to discuss its generalization. 

^ These assumptions can be supported by the perturbative results (13.121) . Also note that as 5 = 0 and h = 0, 
the natural domain of analyticity of the classical dost functions (14.1011 . as functions of the complex variable 
J = is the lower half plane —tt < arg(J) < 0 which corresponds to the strip 0 < 3m(a) < tt. 
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The boundary operators ag are particular representatives of the family of “primary” bound¬ 
ary helds 


^(a) ^ 


(Js 


(s = 0,3,±) 


( 6 . 1 ) 


(here a± = cri±i(T 2 ). Recall that the Hamiltonian ij' (l2.11l) relates to the renormalized Hamilto¬ 
nian H(j, fl2.15l) by means of the canonical transformation HU with U = exp (icr3(/)(0)) 

and $^(0) = 20(0). For this reason the scale dimension of (T±e‘“‘*’® is given by D±{a) = g{a^l)‘^ 
whereas the scale dimension of cro, 3 e^“‘*’® equals to Ho, 3 ( 0 ) = 1?®^- An overall normalization of 
the primary boundary helds can be chosen to satisfy the normalization condition at ?/ —)■ 0^, 

^ + ... (S = 0,3) 

(ii = ±), (6.2) 

along with the hermitian conjugation 






for s = 0, 3 
for s = ± 


(6.3) 


Repeating the construction from Section [T3l one can introduce the normalized partition 
function Z 2 i{p\y) fl2.18p for the operators O = As y ^ O’*", the effect of the boundary 
interaction becomes negligible within the interval (0,?/) and the pair of hermitian conjugated 
operators at the ends of the segment can be substituted by a local insertion dehned by the 
leading term of the operator product expansions fl6.2p . This yields the leading behavior at 
y^0+: 




y Z.ga? j for s = 0, 3 

I (1 — smi) y-‘^9{a-sf _|_ 0 ( 1 )) for s = ± 


(6.4) 


The spin degrees of freedom are freezing out at large distances so that the large-?/ asymptotic 
of Z 21 (cri“^ I?/) has a similar form as eq. (12.231) . 


= l4l2i(a(“))|' y-^d,,{a) (f + o(l)) as ?/ ^ +CX) . (6.5) 

The IR exponent 0 ( 21 ( 0 ) can be found using the mean held approximation from Section lXTl 

0 ( 21 ( 0 ) = ^ {m2 — mi — 2 aY ■ ( 6 . 6 ) 

Just as in the case a = 0, this formula is expected to be an exact result. This supposition is 
supported by a perturbative calculation given in the appendix. As before, A 2 i(crs“^) in fl6.5p can 
be thought as scaling functions depending on the mi, m 2 , and the pair of Kondo temperatures 
HJ", E 2 . In the context of the spin-boson model (11.4p . the hdelity A 2 i(cri“^) represents the 
universal part of the vacuum-vacuum matrix element (H 2 | I “bare” operator 

(II.6p . which requires both UV and IR regularizations. If the integration in (11.41) and (11.61) is 
restricted to the hnite domain L~^ < k < A, then 




(6.7) 
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In the case hi = /i 2 , the calculation of A 2 i[ai°'^) can be reduced to hnding p-vacuum matrix 
elements generalizing (15.1^ : 


{p I Te^{mi) \p) = mi) (52p'-2p,2a+m2+mi • 


( 6 . 8 ) 


The bosonization approach developed in ref. IZH allows one to study the case hi = ^2 —)■ O'*" and 
extend the result from Section[531 It turns out that, for the primary held with 0 < a < 1, 
the hdelity is simply expressed in terms of the function u{a) (I5.30j) : 


.<7(5-a) 


■4apr)L..,.^„.=S+(a) n{a) (ElElY 

where a — log(iI'l/i?f), In the ease of the corresponding fidelities are given by 

= So(«) (6“““ + (-!)• ) (* = 0,3) 


where 


A^°'\a) = u{a) 




dy m(7 + y)m( 7 — y) 


271 coshy m( 7 — a + y)m( 7 — a — y) 


(6.9) 


( 6 . 10 ) 


( 6 . 11 ) 


Note that 5+(a) and iSo(a) in eqs. fl6.9p . flh.lOp stand for some dimensionless functions of a and 
the coupling g, which remained undetermined except their values at a = 0: 


5+(0) = - 


irg 


t ( ^ rn \( ’ o mi 1 


1 2 r( 


2 + 2(l-g). 


( 6 . 12 ) 


Finally, one can consider the hdelities for scaling helds of the form 




where Pn{- ■ •) stands for a differential polynomial of order N built from the boundary held 
dy^B- In this most general setting, the problem is related to hnding the multipoint amplitudes 

{p' \e^"am ■■■T,^{mi)\p) = . . . , mi) 52p'-2p,-2a+m^+-+m^ ■ (6.13) 
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Figure 5: The Anderson-Yuval [9] one-dimensional gas with an external held h and external 
charges ±a at the ends of segment. The arrow shows the direction of the external force exerted 
on the “positive” charges. 


A Appendix: Conformal Perturbation Theory 


Here we consider the large y asymptotic of the normalized partition function .E 21 (ctq^^ ||/) for 
= 0, hi > 0, or equivalently, mi = 1. Note that in this case 

Z2i{a^^^\y) = Z2i(af^||/) , Z2i{a^^^\y) = 0 . (A.l) 

Using the renormalized Hamiltonian fl2.15p . where /i is renormalized coupling corresponding 
to the Kondo temperature E* = E 2 and h = ^ 2 , the ratio Z 2 i{cy^'^\y) = Z 2 i{cy^'’\y)/Zn can 
be represented in the form of the grand canonical partition function of the Anderson-Yuval 
one-dimensional gas of alternating charges on the hnite interval, 

^2i(4“ ||/)|e*=o = / d|/ 2 n-l-- - / d|/i JJ ('^.2) 

' n=0 >^0 >^0 >^0 1=1 

X (^p \ g-2i0(j/2n)g-l-2i0(y2n-l) . . . ^-2\<f>{y2) g2i(/)(?;i) g+2ia0(O) | ^ ^ 


where, in fact, the matrix element does not depend on the choice of p-vacuum, 

(pIe”“"*(!/„)---e”“‘*(!/i)|p> =<5a.+-+..„o '[{(Vi-V■ (A.3) 

i>j 


The parameter h plays the role of the external electric held whereas the exponentials can 

be interpreted as creators of two external charges at the ends of the interval (see Fig.|5]). Then 
one has 


|^(4“^)|L_n = lim fd2i-2ga^^e{i-eo{\)) ( 1 + g„(£) ) . (A.4) 

' U l=}iy^^QQ \ f 


71=0 


Here ^21 is some (a priori unknown) orthogonality exponent, 

eo(A) = —Eo/h > 0 , A = , 


(A.5) 


and 


pi rU2n rU2 ”■ 

qnii) = / du2n / dU2n-l---/ dUi TT 

^0 Jo Jo 

n 

X Y\ \'^2i-U2j-l\~^^ \{u2j -U2i){u2j-l-U2i-l) \ ^ . (A.6) 

2J = 1 


(1 - U2i)u2i-l 
(1 — U2i-l)u2i 
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Consider the first nontrivial coefficient One can show that as £ ^ +cx3, 

gi(£) = r(l -2g)i- r(2 - 2g) (^1 + Aag log (ie^) ) + o(l) , (A.7) 

where 0=1 + log(2) — 'Je —'^(2 — 2g) + 2ag) and ^p{z) = logr(^). Substituting (IA.7P 

in flA.4ll we observe that the term oc is canceled out provided that eo(A) = 1 + 2^^“^ r(l — 
2g) X^ + 0(A'^). Taking into account that m(A) = eo — (1 — g)X i.e., 

m(A) = 1 - 2^5-1 r(2 - 2g) X^ + O(A^) , (A.8) 

one hnds that the cancellation of the divergent term oc A^ log £ requires that 

^21 = ga^ + ga {1 — m) + 0((1 — m)^) , (A.9) 

which is in agreement with fl6.6p specialized for mi = 1, m 2 = m. The above calculation 
suggests that 

|Al 2 i(crS“^)U*^o = A{m,ag) , (A.IO) 

where 

/TTl 

A(m,/) = J—^ (^l + ^a„(/) (1-m)’"), (ATI) 

n=l 

and the hrst expansion coefficient reads explicitly 

ai{l) = -/(! + log(2) - 7 e - ^P{2 - 2g) - ^^(1 + 2/)) . (A.12) 

The function flA.lOp is expected to have a hnite limit as h —)■ O’*": 

£m p2i(+’)U..„ = (£*)-“'<“+!> +“> . (A.13) 

The dimensionless amplitude is the universal part of the partition function of the 
system depicted in Fig.[6l 

~ {cisEy^ (c2LO*)-5(“+i)' . (A.14) 

Here ^free is the partition function of the Gaussian CFT with free BC, e is the lattice spacing, 
and a-independent, non-universal coefficients Cj depend on details of the IR and UV regular¬ 
izations. 

It deserves to be mentioned here that | A 2 i(cro“^) |g*_Q has a hnite limit as g ^ 0 and I = ag 
is kept hxed. In this case, as it follows from flA.2F flA.3F 

Im. y»“^22i(+l!/)|,..„ = (t|r,exp 

/ = ag —fixed 


J dr(i 


,2i0c(-iT) 


(j_ + e 


-2i0c(-ir) 


^+) 


t)(A.15) 
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{E* = 0, h>0) 


Gaussian CFT 


L —)■ oo 




{E* >0, h = 0+) 


Figure 6: The dimensionless amplitude flA.13D is the universal part of the partition function 
of the statistical system schematically shown in this hgure. 


where i0c(—i'r) = hr + I log (^^)- In fhe presence of singularities at r = 0, y, the analysis 
performed in Section lTT] requires some modihcation. One can show that A{m, 1) a.t g = 0 
coincides with the connection coefficient 


T^\m,l) = lini (-r)'Tk(-ir) , 


2 T^0“ 

where Tk(—ir) is the dost solution of the Sturm-Liouville equation 




2k/m /(/ + 1) 


+ 


r 




satisfying the asymptotic condition 


^k(-iT') ^ (-t) 


—Im ^kr 


as 


T = 0 , 


r —oo 


(A.16) 


(A.17) 


(A.18) 


(recall that we use the notations k = \/J'^ + > 0 and m = h/k). The dost solution is 

expressed in terms of the confluent hypergeometric function: 


Tk(-ir) = k'™ (-kr)^+' o’* 
This yields the result 


2 i+;(i+m) ^ ^ ^ 2(1 + 1 ), -2kr) 


/c)/ l + m . 2^(”^-^)Fl + 2 

+ ^ ^ 2 F 1+ l + m 


(A.19) 


(A.20) 


Notice that an expansion of T‘^\m,l) at m = 1 turns out to be consistent with the classical 
limit of flA.10|) - flA.12p . Thus we see that | A 2 i(cro°'^) | b*=o = h~^ and it is an analytic 

9=0 

function for 3?e(/) > —There are two other points which deserve notice, such as the bilinear 


relation 


where 


ri"^(e^"m, 1 ) 1 ) + Tl"^(e‘"m, 1 ) T|"^(m, 0 = 1, 


dc) 


(c), 


dc)/ 




1-^ 2^(-+^)F(l - 21 ) 

2 F(1 - (1 -m)/) ’ 


(A.21) 


(A.22) 






























(to be compared with eq. fl5.10p ). and the behavior of as m —)■ O’*': 

1) m! 2' . (A.23) 
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